On initial condition of the new agegraphic dark energy model 
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The initial condition £2 de (z ini ) = n 2 (\ + Zi n i) _2 /4 at z ini = 2000 widely used to solve the differential equation 
of £2 de , the density of the new agegraphic dark energy (NADE), makes the NADE model be a single-parameter 
dark-energy cosmological model. However, this initial condition, we find, is only applicable in a flat universe 
with only dark energy and pressureless matter. In fact, in order to obtain more information from current obser- 
vational data, such as cosmic microwave background (CMB) and baryon acoustic oscillations (BAO), it often 
needs us to consider the contribution of radiation. For this situation, the initial condition mentioned above be- 
comes invalid. To overcome this shortage, we deeply investigate the evolution of NADE in the matter-dominated 

and radiation-dominated epochs, and obtain a new initial condition £2 de (z ini ) = - (1+ ^ ini) (l + V^foni)) at 
Zini = 2000. Here F(z) = n ~^~^fi) w * m anc * ^mo the current density parameters of radiation and pressure- 
less matter, respectively. This revised initial condition is applicable for the differential equation of Q de obtained 
in the standard FRW universe with dark energy, pressureless matter, radiation and even spatial curvature, and 
can still keep the NADE model being a single-parameter model. With the revised initial condition and the ob- 
servational data of type la supernova (SNIa), CMB and BAO, we finally constrain the NADE model. The results 
show that the single free parameter n of the NADE model can be constrained tightly. 
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I. INTRODUCTION 

The accelerated expansion of current universe, first found in 
1998 [1], implies that our universe is being dominated by an 
exotic component with negative pressure, dubbed dark energy. 
To understand its nature, we should first ascertain its dynamic 
evolution. For many dark energy models, it is believed that 
the models are favored by observations if they can fit the data 
well, and meanwhile they have less free parameters, since the 
less parameters may be constrained tightly. To our knowledge, 
in dark-energy cosmology, there exist three rare dark energy 
models, namely, the Lambda Cold Dark Matter (ACDM), the 
Dvali-Gabadadze-Porrati (DGP) braneworld [2] and the new 
agegraphic dark energy (NADE) [3] models which only con- 
tain one free model parameter. Among these three models, 
the NADE model is a special one, because unlike the for- 
mer two models whose single-parameter feature is obvious, 
the NADE model is due to its special analytic feature in the 
matter-dominated epoch [4] . To see it clearly, we first briefly 
review the NADE model. 

The dark energy density pde of NADE model, constructed in 
light of the Karolyhazy relation [5] and corresponding energy 
fluctuations of space-time, has the form [3] 



Pde = 



3n 2 Ml Y 
rj 2 



(1) 



where n is a numerical parameter, Mpi is the reduced Plank 
mass, and r] is the conformal age of the universe, 
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(2) 



where a is the scale factor of the universe and H = a/ a is 
the Hubble parameter. Here an overdot denotes the derivative 
with respect to the cosmic time t. In a flat universe containing 



dark energy and pressureless matter, the Friedmann equation 
can be written as Dde + = 1 , where Qde and Q m are defined 
as the ratio of the densities of dark energy pa e and matter p m 
to the critical density p cr i t = 3M 2 X H 2 , respectively. From the 
Friedmann equation, Eqs. (1), (2) and the energy conservation 
equation p m + 3Hp m = 0, one can derive a differential equation 
ofQdefe) [3]: 



dz 



'de 



l+z 



d-flde) 



3 V^de 



(3) 



where z = a~ l - 1 is redshift. Furthermore, combining Eqs. (1) 
and (2) with the energy conservation equation pde + 3//(l + 
w)Pde = 0, one can easily find that the equation-of- state pa- 
rameter (EOS) of NADE is given by [3] 



w = -1 + 



2 VHde" 

3n a 



(4) 



At first glance, one might consider that NADE is a two- 
parameter model, since besides the parameter n the model 
has another free parameter Q m o coming from the natural ini- 
tial condition QdeO = 1 _ ^mo in solving Eq. (3) (note that 
the subscript "0" denotes the present value of the correspond- 
ing quantity, hereafter). However, as shown in Ref. [4], the 
NADE model is actually a single-parameter model in prac- 
tice, thanks to its special analytic feature Dde = n 2 a 2 /4 in the 
matter-dominated epoch. To obtain this relation, one can con- 
sider the matter-dominated epoch, in which H 2 oc p m oc a~ 3 . 
From Eqs. (1) and (2), one obtains pde °c a~ l . Then, from the 
energy conservation equation pd e + 3/f (1 + w)pd e = 0, one has 
w = -2/3. Comparing w - -2/3 with Eq. (4), one can find 
that Q de = n 2 a 2 /4 = n 2 (l + z)" 2 /4. Note that this relation is 
also one of the analytic solutions of 



~dz~ 



-a 



de 



l+z 



3 — V^de 



(5) 
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which is the reduced form of Eq. (3), since 1 - Qde - 1 in 
the matter-dominated epoch. Thus, once the value of n is 
given, Eq. (3) can be numerically solved by using £lfe(Zmd - 
n 2 (l + Zini)~ 2 /4 at any zm deep into the matter-dominated 
epoch (zini = 2000 is chosen in Ref. [4]). 

By using the initial condition ^defcni) = n 2 (l + Zini)~ 2 /4 at 
Zini = 2000 and the observational data, Wei and Cai [4] firstly 
constrained the single-parameter NADE model. After their 
work, this initial condition was widely used in the literature; 
see, e.g., Refs. [6-8]. All results showed that the only free 
parameter n could be constrained tightly and the NADE model 
could fit the observational data well. 

However, the fly in the ointment is that the initial condition 
^defcni) = n 2 (\ + Zini)~ 2 /4 at zm = 2000 is obtained by con- 
sidering a flat universe containing only dark energy and pres- 
sureless matter. So, a natural question we may ask is whether 
this initial condition is applicable when we consider the con- 
tribution from radiation. Before we answer this question, let 
us first see why we have the necessity to extend this discus- 
sion. 

We all know that the cosmic microwave background (CMB) 
and large-scale structure (LSS) observations play an essential 
role in testing the cosmological model and constraining its ba- 
sic parameters. Generally, we might need to use the full data 
of CMB (CMB temperature and polarization power spectra) 
and LSS (matter power spectrum) to perform a global fitting. 
However, such a way of fitting consumes a large amount of 
time and power. As an alternative, two methods are widely 
utilized: (i) using the shift parameter R from CMB [9, 10] 
and distance parameter A of the BAO measurement [11], 
and (ii) employing the distance prior including R, l A and z* 
from CMB [12] and r s (zd)/D v (z) from BAO measurement of 
SDSS [13]. Here r s (z) is the comoving sound horizon whose 
fitting formula is given by 



r s (z) 



~ f 

V3 Jo 



da 



a 2 H(a) ^1 + (3Q, b0 /4Q. y0 )a 



(6) 



where the present photon density parameter £2 r o = 2.469 x 
10 _5 /T 2 (for r cmb = 2.725 K) with h the Hubble constant 
Ho in units of 100 km s" 1 Mpc -1 , and Qj,o is me present 
baryon density parameter. For other quantities mentioned in 
the method (ii), we will illustrate them in detail in Sec. III. 
There is no doubt that the method (ii) encodes more informa- 
tion of the CMB and LSS data. Note that the distance prior 
(R, I a, z*), as shown in Ref. [12], is applicable only when 
the model in question is based on the standard FRW universe 
with pressureless matter, radiation, dark energy, as well as 
spatial curvature. Since the integral in Eq. (6) involves the 
early radiation-dominated epoch, it needs us to consider the 
contribution of radiation when utilizing r s (zd)/D v (z) of BAO 
measurement. To sum up, in order to use the method (ii), we 
need to at least consider the contribution of radiation besides 
the pressureless matter and dark energy. 

In the following, we will show that the old initial condition 
^defcni) = ft 2 (l + £ini)~ 2 /4 at zm = 2000 needs to be amended 
to accommodate the method (ii) in using the CMB and BAO 
data. Then, with the revised initial condition and the current 
observational data including type la supernovae (SNIa), CMB 



and BAO, we will constrain the NADE model in a flat universe 
with dark energy, matter as well as radiation. 



II. NEW INITIAL CONDITION 

In a flat universe with dark energy, matter and radiation, the 
Friedmann equation reads 



Qde + Q-m + = 1? 



(7) 



where £l r is the ratio of the energy density of radiation p r to 
the critical density p cr i t . Using Eqs. (1), (2) and (7), and com- 
bining the energy conservation equations p m + 3Hp m = and 
p r + AHp r = 0, one can easily derive the differential equation 

of Qde(z), 



da 



de 



-ft, 



de 



l+z 



d-ftde) 



O j-y \ 2(1 + Z ) 

3 + F(z) Vftje 



(8) 



To solve Eq. (8), of course, we may use the initial condition 
^deo = 1 - ^mo _ ^r0- However, such treatment will add an 
extra parameter D m0 to the NADE model, as mentioned above 
(note that one usually fixes Q r o = £Xyo(l + 0.227 liVeff), and 
the standard value 3.04 of the effective number of neutrino 
species Af e ff is required [12]). Moreover, it has been shown 
in Ref. [8] that such a two-parameter NADE model cannot be 
constrained well by observational data (e.g., 137.702 < n < 
331. 91 A at the l<x level). On the other hand, if we expect that 
^defcni) = ft 2 (l + £ini)~ 2 /4 at Zim = 2000 is able to be used as 
the initial condition in solving Eq. (8), we must require it at 
least to satisfy the equation 



d£lfe 

dz 



-a 



'de 



(9) 



which is the reduced form of Eq. (8). Here 1 - £2de - 1 
in the matter-dominated epoch is used. Comparing Eq. (9) 
with Eq. (5), we can find that the above condition depends 
on F(z) « 1 in the matter-dominated epoch, since £2de = 
n 2 (l + z)~ 2 /4 satisfies Eq. (5) accurately. However, from the 
definition F(z) = Q ^+n*(i +z ) » ^ ma Y nave a visible value be- 
cause of z » 1 in the matter-dominated epoch, even though 
Q r0 «c n m0 . In fact, one can check that F(2000) = 0.3850 if 
we choose Q m o = 0. 1 334 x h~ 2 according to the recent WMAP 
observations [12]. Thus the existence of the non-ignorable 
F(z) indicates that we need to find a new initial condition to 
solve Eq. (8). 

Fortunately, it's not that hard to obtain the new initial condi- 
tion applicable for Eq. (8). Let us consider an epoch when the 
density of dark energy can be ignored; but we don't strictly 
distinguish the matter-dominated or the radiation-dominated 
epoch; and then we have H 2 oc Q m0 ^~ 3 + ^ r o^~ 4 - From 
the definition of the conformal age of the universe, we have 
r] oc VQ m o<z + Ok) - V^rO- Then, from Eq. (1), we can obtain 



Pde OC ( VQ m o<3 + ^rO - V^rO ) 



(10) 
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Combining Eq. (10) with the energy conservation equation 
Pde + 3H(1 + w)pde = 0, we can obtain the EOS of dark energy 
at that epoch, 



2 i a 

3 3 V Qm{)Cl + £2r0 



(ID 



Then, comparing Eq. (11) with Eq. (4), we can easily obtain 



n\\+z)- 



(l + jm) . 



(12) 



Furthermore, one can also check that Eq. (12) is an ana- 
lytic solution of Eq. (9). As the evolution of Dde satisfies 
Eq. (12) both in the matter-dominated epoch and the radiation- 
dominated epoch, we can use Eq. (12) at any zm in these 
two epochs as an initial condition to solve Eq. (8) numer- 
ically. In our work, we follow Ref. [4] and still choose 
= 2000. Then, a new initial condition ^defcni) = 
n2(1+ 4 Zini) " 2 (l + V^XiiZ)) 2 at = 2000 is available for the 
NADE model in a flat universe with dark energy, pressure- 
less matter and radiation. In fact, this new initial condition 
is also valid in a non-flat universe, since the spatial curvature 
£lk is much smaller than Q m or £l r at z = 2000. It is interest- 
ing to make a comparison between the new initial condition 
with the old initial condition Qdefcni) = ^ 2 (1 + £ini)~ 2 /4 at 
= 2000. Obviously, their difference comes from the term 

(l + V^fani)) whose value is 2.626 at zm = 2000, which 
means that the density of dark energy at z = 2000 given by 
the new initial condition is about 2.6 times greater than that 
given by the original initial condition. 

Up to now, we have discussed the new initial condition the- 
oretically. Next step, we might expect, is to test the new initial 
condition observationally, namely, using observational data 
to constrain the NADE model with the new initial condition. 
To achieve this, we will use a Markov Chain Monte Carlo 
(MCMC) method. However, before doing this, we still need to 
overcome a technical difficulty. It's well known that we need 
artificially give a start point of the free parameters to launch 
the MCMC. For our work, as our purpose of trying to find 
the new initial condition is to keep the single-parameter fea- 
ture of the NADE model, we only give an initial value of the 
single free parameter n. Here, the difficulty occurs, which is 
how we can numerically solve Eq. (8) by using the new initial 



condition £2 de (zini) 



» 2 (1+Zini)" 2 



(l + V^XiZ)) at Zini 



2000, 



as we only know the initial value of n but have no idea about 
Q m o (the value of Q r o is fixed as mentioned above). Note that 
both Eq. (8) and the new initial condition explicitly contain 
the parameter £l m o. 

To solve the problem, actually, we have many methods. 
Here, we introduce two methods we use. The first method 
is treating the non-independent parameter £l m o as a variable. 
Thus, Eq. (8) becomes a partial differential equation, 



d In Q d e _ Qde - 1 



dz 



l+z 



3 + 



Q r0 (l+z) 2(1 +z) 
Q m0 + ^ro(l+£) n 




where Qde = ^de(z, ^mo) is a function of the two vari- 
ables {z, Q m o}- Thus, giving a value of n and a range of 
Q m o, one can numerically solve the partial differential equa- 
tion (13) by using the new initial condition Qdefanb ^mo) = 

" 2(1+ 4 Zim) ~ 2 (l + V^fcni, ^mo)f at zm = 2000. Then, with the 
result of £2de(z, ^ m o), we can obtain the value of Q m o by 
numerically finding the root of the equation £2d e (0, Q m0 ) = 
1 - Q m o - n r o- Substituting Q m o back into Qdefe Q m o), we 
can obtain the evolution of £lde(z). 

Besides, we can also use a numerical iterative method, 
namely, generating a sequence {£^ } by the iterative formula 

Q^ 1} = 1 - Qde(0)| n (o - n r0 . Here, n de (z)| a (o is the nu- 
merical solution of Eq. (8) with £l m o = Q® Q . Thus, for the 
current value of oP Q , the new value can be obtained 

from the iterative formula. We can choose the initial value 



Q 



(0) 
m0 



0.27 and set a termination condition for the iteration, 



such as IQ^q 1 ^ ~~ ^mol - 6 w ^ m 6 a small quantity. Our prac- 
tice shows that its convergence speed is very fast, and gener- 
ally, the number of the iteration is less than 5 for e = 10" 6 . 

Using the two methods mentioned above, we finally find 
the evolution of Qdefe) for a given value of the single free 
parameter n. Then, the dimensionless Hubble expansion rate 
is given by 



E(Z) = 



H(z) 
Ho 



D m0 (l+z) 3 +n r0 (l+z) 4 



1 - n de (z) 



1/2 



(14) 



III. OBSERVATIONAL DATA AND RESULTS 

In this section, we constrain the NADE model with the new 
initial condition by using the data from Union2 SNIa (557 
data) [14] and observations of CMB from 7-year WMAP [12] 
and BAO from SDSS DR7 [13]. 

The data points of the 557 Union2 SNIa are compiled in 
Ref. [14]. The theoretical distance modulus is defined as 



Ah(z/) = 5 log 10 D L (zd + yU , 



(15) 



where //o = 42.38 - 5 log 10 h and the Hubble-free luminosity 
distance 



Z>l(z) = (1+z) 



r dz' 
Jo E(z\ 



oy 



(16) 



where 6 denotes the model parameters. Correspondingly, the 
X 1 function for the 557 Union2 SNIa data is given by 



557 2 

x 2 (6) = ^ !^ obs( ^ 



(17) 



where cr is the corresponding lcr error of distance modulus 
for each supernova. The parameter //o is a nuisance parameter 
and one can expand Eq. (17) as 



X 2 SN (0)=A(0)-2fi o B(0) + fi 2 o Q 



(18) 
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where A(0), B(0) and C are defined in Ref. [15]. Evidently, 
Eq. (18) has a minimum for /jq = B/C at 



X 2 SN (0)=A(0)- 



B(0) 2 
C 



(19) 



We use the BAO data from SDSS DR7 [13]. The distance 
ratio (d z ) at z = 0.2 and z = 0.35 are 



do.2 



r s (Zd) 
D V (02Y 



35 



r s (Zd) 
D V (035Y 



(26) 



Since ^ N min = X SN ,min> instead minimizing^ we will min- 
imized^ which is independent of the nuisance parameter j^o. 

For the observational data of CMB and BAO, we have men- 
tioned in Sec. I that two simple but efficient methods are often 
used instead of using their full data to perform a global fit- 
ting. Among these, the method (ii) using distance prior (R, 
l A , z*) from CMB and r s (z d )/D v (z) at z = 0.2 and 0.35 from 
BAO, contains more information of CMB and BAO observa- 
tions but requires considering the contribution from radiation. 
Since our new initial condition is mainly designed for this re- 
quirement, we will adopt the method (ii) in our work. 

The "WMAP distance prior" is given by the 7-year WMAP 
observations [12]. This includes the "acoustic scale" l A , the 
"shift parameter" R, and the redshift of the decoupling epoch 
of photons z*. The acoustic scale l A describes the distance 
ratio D A (z*)/r s (z*), defined as 



l A = (1 + z *) 



(20) 



where a factor of (1+z*) arises because D A (z*) is the proper an- 
gular diameter distance, whereas r s (z*) is the comoving sound 
horizon at z* and its fitting formula is given by Eq. (6). For 
the value of Qj,o, we fix = 0.02246 x /T 2 , given by the 
7-year WMAP observations [12]. We use the fitting function 
of z* proposed by Hu and Sugiyama [16] 

z , = 1048[1 + 0.00124(Q^ 2 )" a738 ][l + gii^moh 2 ) 8 2 ], (21) 



where 



2\-0.238 



0.0783(Qm)/z 2 ) 
1 + 39.5(£W* 2 ) 



'.763 ' 



§2 = 



0.560 



1 + 2l.l(Q b0 h 2 ) lM ' 

(22) 

The shift parameter R is responsible for the distance ratio 
D A (z*)/H-\z>), given by [17] 

R(z*) = ^moH 2 (l + z*)D A (z*). (23) 

Following Ref. [12], we use the prescription for using the 
WMAP distance prior. Thus, the x 2 function for the CMB 
data is 



(24) 



where xi - (l A ,R,z*) is a vector, and (C^^ij is the inverse 
covariance matrix. The 7-year WMAP observations [12] give 
the maximum likelihood values: l A (z*) = 302.09, R(z*) = 
1.725, and z* = 1091.3. The inverse covariance matrix is also 
given in Ref. [12]: 



( 2.305 29.698 -1.333 ^ 
29.698 6825.27 -113.180 
-1.333 -113.180 3.414 



(25) 



where r s (zd) is the comoving sound horizon at the baryon drag 
epoch [18], and 



1/3 



(27) 



encodes the visual distortion of a spherical object due to the 
non Euclidianity of a FRW spacetime. The inverse covariance 
matrix of BAO is 



i _ / 30124 -17227 
^bao) ~ | -17227 86977 



The^ 2 function of the BAO data is constructed as: 



-1 



jobs\ 



(28) 



(29) 



where d\ - (do.2, do.35) is a vector, and the BAO data we use 
are d 02 = 0.1905 and d 03 5 = 0.1097. 
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FIG. 1: The probability contours at lcr and 2<x confidence levels in 
the n — h plane for the NADE model with the new initial condition. 



The best-fitted parameters are obtained by minimizing the 
sum 




2 _ ~2 2 2 

X - /tsN + ACMB + XBKO ' 



(30) 



Using the MCMC method, we finally find the best-fit param- 
eters: n = 2.810^ }}3 an d h = 0.665+g-gJ| at the lcr level, 
and n = 2.810+g;}j*'and h = 0.665+°;°^ at the 2cr level. The 
best fit gives values of^^ = 577.451 and Q m o = 0.268. In 
Fig. 1, we plot the contours of lcr and 2<x confidence levels 
in the n - h plane for the NADE model with the new initial 
condition. 

The obtained results show that with the new initial condi- 
tion, the observational constraints on the single parameter n 
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of NADE model are fairly tight. Next, we would like to com- 
pare the results of this work with those of previous works. We 
choose the work in Ref. [8] to make a comparison. Note that 
the work in Ref. [8] is also about the observational constraint 
on the NADE model but based on the old initial condition, 
and the observational data used in Ref. [8] also come from 
the Union2 SNIa, 7-year WMAP and SDSS DR7 except that 
the method of using CMB and BAO observations is different 
from our present work (method (i) in Ref. [8] while method 
(ii) in our present work). So, the fitting results in these two 
works are comparable. The best-fitted n and the correspond- 
ing Q m0 in that work are 2.886 and 0.265, respectively. We 
see that though the values of n in the two works are relatively 
shifted by a small number, the values of Q m o produced by 
the model are rather similar. For the 95% limits on n, the 
work of Ref. [8] gives 2.723 < n < 3.055, and the present 
work gives 2.625 < n < 2.996. To see the influence of the 
values of n in the two works on the model, we will make a 
comparison with the quantity f(z) = £lfe(z)l£l m (z). Using 
the best-fit results obtained in the two works, we find that 
/(0)new//(0) o id = 0.986 and /(2000) new //(2000) old = 4.362. 
So, it is clearly seen that though the difference in the early 
epoch is fairly evident, the results produced in today are sim- 
ilar. 



IV. CONCLUSIONS 

The NADE model is considered to be a single-parameter 
model, due to its special analytical feature Qde = n 2 (l +z)~ 2 /4 
in the matter-dominated epoch. Thus, once the value of the 
single free parameter n is given, the differential equation of 
£lde(z) can be numerically solved by using the initial condition 
^defcni) = ft 2 (l +£ini)~ 2 /4 at Zmi = 2000. However, this initial 
condition is only applicable in a flat universe with only dark 
energy and pressureless matter. That is to say, when we need 
to consider the contribution from radiation, this initial condi- 
tion becomes invalid. On the other hand, some cases indeed 
need us to consider the contribution of radiation. For instance, 
when using the CMB and BAO data to constrain dark energy 
models, we mainly have two methods — method (i): using the 
shift parameter R from CMB and distance parameter A from 
BAO, and method (ii): employing the distance prior includ- 
ing R, l A and z* from CMB and r s (zd)/D v (z) from BAO. Of 



course, the method (ii) encodes more information of the CMB 
and BAO observations. However, the method (ii) requires us 
to consider the contribution from radiation. 

Thus, in order to utilize the method (ii) to fit CMB and BAO 
data, we deeply analyzed the NADE model in a flat universe 
with dark energy, matter and radiation. Finally, we found a 

similar analytical solution, Dde = n (1 ^ z) (l + V^Xi)) with 
F(z) = Q + (i + z) > m the early epoch (matter-dominated 
or radiation-dominated epoch). For the initial z, we still 
chose Zini = 2000. Hence, we have a new initial condition 

^defoni) = " 2(1+ 4 mi) ~ 2 (l + ^^F(^j) 2 at zm = 2000. Further- 
more, we found that this new initial condition is also applica- 
ble in a non-flat universe. 

For the problem how we can solve the differential equa- 
tion of Qde before we know the value of the non-independent 
parameter Q m o, we provided two methods to overcome this 
technical difficulty. The first method is to consider the non- 
independent parameter Q m o to be a variable and use the new 
initial condition to solve the partial differential equation of 
Qde- The second method is a numerical iteration method. 
Our practice shows that its convergence speed is very fast. 
With the two methods, we have constrained the NADE model 
by using the new initial condition and the observational data 
including the Union2 SNIa, CMB from 7-year WMAP and 
BAO from SDSS DR7. Our fitting results show that the val- 
ues of n and h can be constrained tightly: n = 2.8 10^- J}^, 
h = 0.665+°;°}* at the la level, and n = 2.810 + °;}^, h = 
0.665+q°2j at the 2<x level. The best fit gives values of 
y 2 . = 577.451 and Q m0 = 0.268. We believe that our new 
initial condition will play a crucial rule in constraining the 
NADE model in the future works. 
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